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Abstract 

A method to construct fractal surfaces by recurrent fractal curves is provided. First we 
construct fractal interpolation curves using a recurrent iterated functions system(RIFS) with 
function scaling factors and estimate their box-counting dimension. Then we present a method 
of construction of wider class of fractal surfaces by fractal curves and Lipschitz functions and 
calculate the box-counting dimension of the constructed surfaces. Finally, we combine both 
methods to have more flexible constructions of fractal surfaces. 
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1 Introduction 

A Fractal surface (or fractal curve) is a fractal set which is a graph of some continuous function on 
(or R). The method of construction of fractal surfaces (or fractal curves) are closely related to 
the generation of fractal interpolation functions(FIF). The FIFs were introduced by Barnsley [3] in 
1986 and after that have been widely studied and used in approximation theory, image compression, 
computer graphics and modeling of natural surfaces such as rocks, metals, planets, terrains and so 

on [II 01 unmans]. 

The constructions of fractal surfaces, such as self-similar, self-affine or non self-affine surfaces, 
by IFSs or RIFSs have been studied in many papers (see [TH HU |SJ [THl HH H]). These surfaces 
are all attractors of some IFSs or RIFSs. In particular, Malysz [T^], Metzler et al [15] and Feng 
et al [S] studied the construction of fractal interpolation surfaces on arbitrary data sets using IFS. 
Furthermore, Malysz [T^], Metzler et al [T5] and Yun [TH] estimated the dimensions of the result 
surfaces. In [5], the authors studied the construction of recurrent fractal interpolation surfaces. In 
[3 [12] they use constant contraction factor in construction of IFS or RIFS and in [15j [HI [9] they use 
function contraction factor in construction of IFS. 

The fractal properties of such rough surfaces as those of metals or rocks may expressed by sectional 
profiles of those surfaces. In many papers, constructions of fractal surfaces by fractal curves have been 
studied. In [T^, Mandelbrot suggested that the fractal dimension of a surfaces constructed by a single 
curve can be obtained by adding 1 to the fractal dimension of the curve. In ISj, Falconer introduced 
fractal surfaces constructed by the movement of a fractal curve along a segment and determined 
their fractal dimension. Xie and Fang proposed the so called star product fractal surfaces which are 
constructed by the movement of a fractal curves along another one(P3)- ^ construction of 



fractal surfaces by 4 fractal curves, which are boundary curves of the constructed surface, and in [TB] 
a construction of Bush type fractal surfaces by two Bush curves were studied. In ff] a construction of 
fractal surfaces by the fractal interpolation. All the constructions of fractal surfaces by fractal curves 
have common property that the fractal dimension of the constructed surfaces is determined by one 
of the fractal curves constructing them. 

We construct recurrent fractal curves by more general RIFSs with function scaling factors and 
estimate the box-counting dimension of the constructed curves. And we construct wider class of 
fractal surfaces by fractal curves and Lipschitz functions, and calculate the box-counting dimension 
of the constructed surfaces. Finally we combine both constructions to construct fractal surfaces. 



2 Construction of Recurrent Fractal Interpolation Curves 

In this section, we construct recurrent fractal curves and estimate their box-counting dimension. 

Let a data set be P = {{xi,yi) e = 1, . . . ,n] , {xq < xi < . . . < Xn) and let Nn — 
{1, . . . , n} , / = [xq, Xn],Ii — [xi-i.Xi],i — 1, . . . , n. We denote Lipschitz(or contraction) constant of 
Lipschitz(or contraction) mapping / by (or cj). Let Z > 2,^ S N and let = [^^(fc), a;e(/c)], a;s(/j), 
Xe(k) S {a;o, . . . , Xn}, here e{k) - s(fc) > 2, fc = 1, . . . , Z . 

For each i E Nn , we fix a fc S {1, . . . , and let k — ^{i) . For i G Nn , k = 7(1) , let a mapping 
Li^k '■ Ik ^ li be contraction homeomorphism and satisfy L^jt : {^Xs(k)TXe{k)} {xi-i,Xi\, and a 
function Fi^k Ik ^ R ^ R defined by Fi^}.{x,y) — Si^k{Li,k{x))a{y) -\- b{x) and satisfy 

Fi^k{xa,ya) ^Vii.aE {s(fc),e(fc)} , (1) 

where Li^k{xa) = xp^jS G {i — and a{x),b{x) are Lipschitz mappings and Si_k(x) is a contraction 
function on li with \si^k{x)La\ < 1 (which is called a vertical scaling factor.) An example of Fi^h 
satisfying (1) is 

Fim{x, y) = Si^k{Li.k{x)){a{y) - g{x)) + h{Li^k{x)), 
where h{x),g{x) are Lipschitz mappings on / and satisfy the following conditions; 

9{xa) ^ya,a& {s{k),e{k)} , h{xi) = y.,J = 0, 1, . . . ,n, 

and Si^kix) are taken as free unknown function. 

We define transformations Wi^k : /fe x R ^ x R (i = 1, . . . , n; k — ^{i)) by 

Wi^k(x, y) = {Li^k{x), Fi^k{x, y)). 

Then there exists some distance equivalent to the Euclidean metric on R2 such that Wi^k{i = 
k — 7(i)) are contraction transformations with respect to the distance. This family 
{Wj^fc : i — 1, . . . ,n} oi transformations defines a recurrent iterated function system([4^) on R^ 
and it has the unique attractor. 

We define a row- stochastic matrix M = {pij)nxn by 




0, otherwise 



where the number indicates for every fixed i, how many = l,...,n include li. Then a 

connection matrix C — {cij)nxn is defined as follows 

r 1, PJ^ > 0, 

1 0, = 0. 
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It is clear that if the row-stochastic matrix is irreducible, then the connection matrix is also irre- 
ducible. 

We denote the attractor of the recurrent iterated functions system(RIFS) {R^;M, Wi,fc,i = 
l,...,n, 7(1) = k G {!,...,/}} by A. Then the following theorem shows that >l is a recurrent 
fractal curve. 

Theorem 1 The attractor A constructed above is a graph of some continuous function which inter- 
polates the data set P. 

(Proof) Let C(J) = {(p G C^{I); (p{xi) = yi,i = 0, 1, . . . , n}, then the set C(/) is complete metric 
space with respect to norm || • ||oo- We can see easily that the operator T : C(/) — )■ C(/) ; (T(p)(x) = 
Fj,/;(ir^(a;),(^(L7fc(^)))'^ e li is well defined and the operator T is a contraction on the complete 
metric space C(/). Therefore the operator T has a unique fixed point in C(/), which we denote by 
/. Then the / is presented by 

f(x) = s^^k{x)f{LJ^{x) + b{x),i • ,n,7(i) = k 

which means that Gr{f) — A , where Gr{f) denote a graph of /. □ 

We estimate box-counting dimensions of the recurrent fractal curves constructed above. We can 
assume that / = [0, 1], since the box-counting dimension is invariant under bi-Lipschitz mapping. 

For a set D{c or R-^) and a function / defined on D, 

Rf[D] = sup{\f{x2)- f{xi)\;xi,X2 € D} . 

is called the maximum variation of f on D. 

Let 7 be a interval in R , L : / ^ / a contraction homeomorphism, a, 6 : R R Lipschitz 
mappings and s : / ^ R contraction mappings with |s(a;)-La| < 1 . We define a mapping F : I xR ^ 
Rby 

F{x, y) = s{L{x))a{y) + b{x), {x, y)GlxR. 

Lemma 1 Let f : I ^R be continuous function. Then 

Rf(l-\ foL-^)[L{I)] < sLaRf[I] + \I\{csaf + Lb). 
Here \I\ is a length of the interval I, s = max/ \s{x)\ and af = max/ |a(/(a;))| . 
(Proof) For any x{€ L{I)), let denote x = L-'^{x){€ I). Then 

\F{L-\f o L-'){x) - F{L-\f o L-^){x')\ = 

= \F{L-\x), f o L-\x)) F{L-\x'),f o L-\x'))\ 
= \s{x)a{f{S:)) + h{i) - s{x')a{f{i')) - b{i')\ 

= \s{x)a{f{x)) - s{x)a{f{x')) + s{x)a{f{x')) - s{x')a{f{x')) + b{x) - b{x')\ 

< sLaRf[I] + Cs\x - x'\af + Lh\x - x'\ 

< sLaRf[I] + \I\{csaf + Lb). 

□ 

Let x^+l - X, = l{i = QA,...,n- 1), .Xe(fe) - x,(k) = ^(A; = 1, . . . J; a e N) and U^k{i e Nn, k = 
j{i)) similitude contraction transformations. Then the number of Ij contained in is a . Let S and 
S_ be diagonal matrices 

S = diag{si, . . . ,s„), S = diag{s^, . . . ,s„), 

respectively, where Si = max/, js, ,y(j)(a;)|, = min/. |.Sj^(j)(a;)|. We assume that the row-stochastic 
matrix M is irreducible and the mapping a(x) is identity. 
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Theorem 2 Let A be the recurrent fractal curve in Theorem [7| // there is some interval Ikg such 
that the points of PD {Ik„ x R) are not collinear, then the box-counting dimension divciBA of A has 
the following lower and upper hounds; 
1) Ifh < 1, then 

1 + loga A < diiiiB^ < 1 + log^ A , 



2)IfX< 1, then 



diiiiB^ = 1 , 



where X — p{SC) and X = p{SC) are respectively spectral radii of the irreducible matrices SC and 
SC. 

(Proof) Proof of 1): Let / be a contraction function whose graph is A. We denote Rf[Ik] by Rk 
and ^ by for simphcity. Then r — )■ oo <^=^ — > . 

After applying each Wi^k to the interpolation points in the interval Ik one time, we have (a + 1) new 
points in every interval Ii{i £ iV„). According to the hypothesis, the interpolation points lying inside 
Jfcg are not collinear and the connection matrix C is irreducible. Thus for every interval li there are 
three points which are not collinear, the maximum vertical distance (computed only along the y-axis) 
from one of the three points to the line through other two interpolation points is greater than 0. The 
maximum value is called a height and denoted by Hi . 
By Lemma [T] on each region li we have 

Rf[Ii\ <s^Rk + ^e, 

where e = Cgf + Li,. 

We define non negative vectors hi, r, Ui and i by 



hi = 



" Hi ' 








" 1 " 


Hn 








1 



Ul 



-ei. 



Since A is the graph of a continuous function defined on /, the smallest number of e^-mesh squares 
necessary to cover x R n ^ is greater than the smallest number of £j.-mesh squares necessary to 
cover vertical line with the length Hi and less than the smallest number of Er-mesh squares necessary 
to cover a rectangle li x [/., fi], where 



/, = min|/(x, y)\, f, = max|/(a;, y)\. 



Therefore, 



i=l 



n 

[H,e;^] - 71 < N{er) < 

i=l 



-ei 



n 



where - > and 



$(hie,Ti) - n < N{er) < $(uie-i + i) 
$(a) = ai + • • • + a„, a = (ai, . 



,a„). 
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After applying Wi^k twice, we have a subintervals of length ^ in each € Nn). Those subin- 
tervals are mapped by the transformation Wj^fe from subintervals lying inside the intervals /^(j) 
corresponding to the interval li containing themslves, and thus for each subinterval /j the height on 
the new subintervals produced in is not less than ■ H, where H is the height on original subin- 
terval contained in the interval Therefore, the sum of maximum variance of / on a subintervals 
of the length ^ contained in the interval G N^) is not greater than i-th coordinate of vector 
U2 = SCui + ^ei, the sum of the heights is not less than i-th coordinate of vector h2 = SChi and 



*(h2£;:^) - an < N{er) < ^{u2e~^ + ai) ^ 



an 



where ^ > £r. 

By induction after taking k such that aSr > ai'-^n — ^^at is, r — log^^ n + l>k>r — log^ n 
and applying k times, we get a''~^ subintervals of the length ^k-in contained in each interval 

and 



^{hks;^) - a'^-^n < N{er) < ^{ukS-^ + a'^-H) 



+ 1 , 



(2) 



where Ufe = SCuk-i + ^bi, = SChfc_i. Then we have 



Ufe = (SC)''-^r + (5C)'=-i -ei + (SC)''-^-ei + ... + (SC)-ei + -ei, 

n n n n 

hfc = (5C)('=-i)hi. 

Since SC and SC are non-negative irreducible matrix, from Probenius' theorem there exist strictly 

positive eigenvectors of SC and 5C(which correspond to eigenvalues A = p{SC) and A = p(SC) 
of SC and SC) and we can choose such strictly positive eigenvectors e , e (which correspond to 
eigenvalues A , A respectively) that 



Then 



r < e, 6i < ne, < < hi . 



N{sr) < ^{uhs;^ + a'=-4) 
<*(ufee-^+a'=-M)(a + l) 

< ^iiSC)''-^rs-^ + {SC)''-^-eieZ^ + (SC)''-^-eis-^ + ... + (SC)-eieZ^ + 

n n n 

n 

< ^{{SCf-^ee;^ + {SC)''~^ee^^ + {SC)''-^ee;^ + ... + {SC)es;^ + 
+ ee7^ +a'="4)(a + 1) 

= (A'=-^$(e)e-i + A'=-^$(e)£-^ + A'=-2$(e)£-^ + . . . + A$(e)£-i + 

< (A'-"$(e)£;i + A'-"$(e)£;i + y-'-^^e)e;^ + ... + X^e)s;^ + 
+ $(e)£7i+a'^-"*(i)(a+l), 

where := log^n. 

On the other hands, since (Sj < {si for i = 1, . . . ,n , from Frobenius's theorem we have A < A. 
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Therefore if A > 1 , then 1 > ^ > i , we obtain 



N{er) < A'-'^$(e)e;i (^l + l + ^ + . . . + ^ + ^ ^^—-] {a + 1) 

= A'-e7iA-''$(e) 1 + + ^ (a + 1). 

Here let S{r) := A-'^$(e) M + ^^^jK + A^-,4(e)a'J + 1)' t^^^n S{r) > and 



By (|2|, we have 



diuisA - hm ^^M^^ < 1 + log^A. (3) 

£r-*-0 —logEr 



Since A > 1 , there is tq such that ri{r) := A^''^^"I>(e) — > for any r(> tq) and therefore we 
have 

^^^^>l + bg,A+-bg,77M, r{>ro) (4) 
By (|3|, Q if A > 1 , then we have 

1 + fog„A < dimsA < 1 + fog<jA. 

Proof of 2): If A < 1, we have 

N{er) < (A'-''$(e)£,7i + y-''<i>{e)e~^ + V~''~''S>{e)e-' + ... + X<^>{e)e^^ + 
+ $(e)£;i+a'-''n)(a + l) 
< e7^[$(e)(r- I/ + 2) +a"''n](a + 1). 

Hence, we have 

diniBA = hm , ^ < 1 + -fog„$(e)(r - ly + 2) + a"''n)(a + 1). 

On the other hands is a curve in and therefore dims^ > 1. Hence dims^ = 1 . □ 

Remark 1. In the case that Si^ k{x) — Si_ ^(constant), ifA = A = A>l, then diniBA = 1 + foga^- 
This is the estimation of box-counting dimension of RFISs in [4] . 

3 Constructions of Fractal Surfaces using Recurrent Fractal 
Interpolation Curves 

In this section, we first introduce some constructions of fractal surfaces in which general fractal curves 
are used and compute their Box-counting dimension. Next, this construction and recurrent fractal 
interpolation curves are combined. 



6 



3.1 Fractal Surfaces Constructed by General fractal Curves and their Box- 
counting Dimension 

We consider fractal surfaces only on the unit square [0, 1] x [0, 1] because our results can easily be 
generalized to any square [a, b] x [c, d]. Let denote / = [0, 1] and E = / x J. 

Let /, g : / — R be fractal curves (i.e. / and g are continuous and their graphs are fractal sets). 
Let A, /i : / — > R be continuous Lipschitz functions with Lipschitz constants L>^ and L^ respectively. 
For {x, y) G E, we define a continuous function : E — > R by 

F{x,y)^\{x)f{x) + f,{y)g{y). (5) 

The following theorem shows that the graph of this function F is a fractal set. 

Theorem 3 Let the function F be given by Then the Box- counting dimension diniB Gr (F) of 
its graph is 

diniBGr (F) = 1 + Max {dimeGr (/) , dimeGr {g)} . 

To prove the theorem, we need some basic concepts and lemmas. Divide the interval / into n 
subintervals with the same length, denote the i-ih subinterval by — [^^, ^] and denote Ejj — 

li X Ij. 

In general, the box-counting dimension dimB.4 of a fractal set A is defined by 

logjV,(^) 

dmiB^l = nm , — 

(5^0 — log 

(if this limit exists), where Ns{A) is any of the foUowings (see [B]): 

(i) the smallest number of closed balls of radius S that cover A; 

(ii) the smallest number of cubes of side S that cover A; 

(iii) the number of (5-mesh cubes that intersect A] 

(iv) the smallest number of sets of diameter at most S that cover A; 

(v) the largest number of disjoint balls of radius 6 with centers in A. 

Let denote Ns {A) on a subinterval Ii by Ng {A) and on a subdomain E^ by Ng-' (A) . In the 
calculation of the Box-counting dimension, we use ^-mesh cubes for n G N {n > 2)(case (iii)). 
The following lemma is easily proved from the definition of the Box-counting dimension. 

Lemma 2 Let A and B be fractal sets and dimB.4 > cVim^B. Then 

^.^log[N.(^)+N.(^)]^ 

5^0 -logo 

log[N^ (^) + N, {B)±5-^] 

lim ; — dmiB^l, 

(5-j-O — log 6 

y log[N^ [A) ±5-^] 

iim ; ; = dimB.4. 

<5-i-0 - log 

Lemma 3 Let the functions /, g be fractal curves. If we define a function F' : E — > R fty 

F'{x, y)^ f{x)+g{y),{x, y) € E, 
then the Box-counting dimension dimBGr {F') is as follows: 

diniBGr {F') = 1 + Max {dimeGr (/) , dimeGr [g)} . 
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(Proof) Evidently Kp' [^ij] — R-/ Pi] + Rg [Ij] for any i, j . On the other hand, it is obvious that 

R^- [E,,] / (^i) < N| {Gr [F')) < 2 + R^, [E,,] / (^i) , (6) 

for the continuous function F' on the domain E^j . Therefore 

Ni (Gr (/)) + Ni (Gr (g)) - 4 < N| (Gr (^^')) < (Gr (/)) + Ni (Gr (.g)) + 2, 
nN 1 (Gr (/)) + nN i (Gr (g)) - in^ < N i (Gr {F')) < 7iNi (Gr (/)) + nN i {Gr (g)) + 2n^. 
From Lemma [2] and the definition of the Box-counting dimension we get the result. □ 

Let denote Mf — Max^gij/ (x) |, cx^ = Miuj-gij A (x) | and c^^ — Max^gijA {x) \ . 

Lemma 4 Let the functions /, A be the same as the above ones. If we define the function Xf : / ^ R 
by (A/) (x) = X{x) f (x) ,x ^ I, then the Box-counting dimension dim^Gr {\f) is the same as that 
of the function f . 

(Proof) For any x, x' S /i, 

\{\f){x)-{\f){x')\ = \\{x)f{x)-\{x')f{x')\ 

= |A (x) f{x)~X (x) f (x') + A (x) / (x') - A (x') / (x') | 
= |A(x)(/(x)-/(x')) + (A(x)-A(x'))/(x')|. 

We can suppose c\-^ > (see Remark 2). Then, for n large enough, 

|A (x) 11/ (x) - / (x') I - La^I/ (x') I < I (A/) (x) - (A/) (x') I < |A (x) ||/ (x) - / (x') I + L,i|/ (x') |, 



cx.Rf [h] - La^M/ < Ra/ \h\ < cx.Rf [/,] + La^M/. 



Thus 



cx,Nl{Gr{f))-m, < Nl{Gr{\f))<cx,Nl[Gr{f))+m2, 
cx,Ni_{Grif))-n-m, < iVj. (Gr (A/)) < ca.Nj. (Gr (/)) + n • m2, 

where mi — 2c\^ + L\Mf, 7712 = 2 + L\Mf. Taking logarithms, the definition of the Box-counting 
dimension and Lemma 1 give the result. □ 

Remark 2 In the case where c\-^ = 0, we choose A' (x) = A (x) -I- c such that cy > and define a 
function A'/ : I ^ R by (A'/) (x) = (A/) {x) + cf (x) for x e L Then (A/) (x) = (A'/) (x)-F(-c/ (x)) . 
Therefore, from Lemma[2]and Lemmajsj we have dimeGr (A/) — dimeGr (/). 

(Proof of Theorem 3) Lemma [s] and Lemma |4] give the result of the theorem. □ 

Corollary 1 Let fi, gj be fractal curves and Xi, fij be Lipschitz functions for i — 1, . . . , A^, f — 
l,...,Af. Then the Box-counting dimension dime Gr (i^) of the function : E — > R defined for 
(x, y)^^ by 

F (x, y) = E^iA, (x) /, (x) -f S^^li/x, (y) g, {y) 

is as follows: 

dimeGr (F) = 1 + Max,,, {dimeGr (/,) , dinieGr {g^)} . 
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The following theorem is proved analoguesly to Theorem [3j 

Theorem 4 Let the functions fi, gj be the same as in Theorem^ and functions Xi, /ij : E — > R 
continuous Lipschitz functions for i — I, . . . , N, j — 1, . . . , M . Then the function : E — > R defined 
for {x,y) GE by 

F {x, y) = SiliA, {x, y) /, [x) + {x, y) g, (y) (7) 

has the Box- counting dimension of 

dimeGr (F) = 1 + Max,,, {dimeGr (/,) , dimgCr (g,)} . 

Remark 3 The fractal surfaces presented in the papers [H [121 113 HO] are contained in the family of 
fractal surfaces defined by ([5]), ([t]). 

3.2 Construction of Fractal Surfaces by Recurrent Fractal Interpolation 
Curves 

The results of the calculations of the Box-counting dimension in Section 3.1 show that the complexity 
of the fractal surfaces defined above are dominated by the fractal curves generating them. Thus, the 
more flexible a construction of fractal curves is, the more natural the fractal surface constructed by 
them is. We can control the complexity and shape of the fractal surfaces in Section 3.1 by the vertical 
contractive function s (x), the stochastic matrix P and Lipschitz functions g{x), h{x), \{x) and ^{x) 
as we want. So, the fractal surfaces presented in this paper should be more appropriate to model 
natural objects. 
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